^ ; The replica symmetric region in the 

(N : Sherrington-Kirkpatrick mean field spin glass 

9^: model. The Almeida-Thouless line 

00 . 

CN ■ Francesco Guerra* 

Dipartimento di Fisica, Universita di Roma "La Sapienza" 

P3 , INFN, Sezione di Romal, Piazzale A. Moro 2, 00185 Roma, Raly 



c/3 



cd 



O 
O 



X 



February 6, 2008 



'^ ■ Abstract 

C ; 

^ , In previous work, we have developed a simple method to study the 

behavior of the Sherrington-Kirkpatrick mean field spin glass model 
for high temperatures, or equivalently for high external fields. The 
^ . basic idea was to couple two different replicas with a quadratic term, 

■^ \ trying to push out the two replica overlap from its replica symmet- 

ric value. In the case of zero external field, our results reproduced 
'!;:j- . the well known validity of the annealed approximation, up to the 

^^ \ known critical value for the temperature. In the case of nontrivial 

f^ ' external field, our method could prove the validity of the Sherrington- 

Kirkpatrick replica symmetric solution up to a line, which fell short of 
2 \ the Almeida-Thouless line, associated to the onset of the spontaneous 

replica symmetry breaking, in the Parisi Ansatz. Here, we make a 
'^ ■ strategic improvement of the method, by modifying the flow equa- 

C , tions, with respect to the parameters of the model. We exploit also 

previous results on the overlap fluctuations in the replica symmetric 
region. As a result, we give a simple proof that replica symmetry 
holds up to the critical Almeida-Thouless line, as expected on phys- 
ical grounds. Our results are compared with the characterization of 
C^ ' the replica symmetry breaking line previously given by Talagrand. We 

outline also a possible extension of our methods to the broken replica 
symmetry region. 
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1 Introduction 

It is very well known, on physical grounds, that for the mean field spin 
glass model, introduced by Sherrington and Kirkpatrick in P, the replica 
symmetric solution holds in a region of the parameters characterized by high 
temperature and/or high external magnetic field, as shown for example in |2] , 
and references quoted there. Recently, Michel Talagrand [3] has been able 
to give a complete characterization of this region, through a method which 
exploits the broken replica symmetry bounds, established in |4j, extended 
to the case of two coupled replicas. Talagrand method is remarkable, not 
only because it gives a full characterization of the replica symmetric region, 
but also because it can be easily extended, with some technical effort, to the 
broken replica symmetry region. In this way, it has been possible to establish 
the rigorous validity of the Parisi Ansatz for the free energy density in the 
infinite volume limit [3|. 

The motivation of this paper can be easily explained. It is surely true that 
some features of the broken replica symmetry phenomenon can be found also 
in the region of the parameters where replica symmetry still holds. This is a 
general aspect of statistical physics. Let us consider for example a system in 
its liquid phase near to the onset of crystallization. While the liquid phase is 
translation and rotation invariant, one can see features of the crystal phase, 
where translation and rotation invariance are broken, through a study of 
the correlations, near the transition line. In the approach of Talagrand, the 
broken replica symmetry features, for coupled replicas, are exploited in an 
essential way in order to characterize the boundaries of the replica symmetric 
region. However, it should be possible, in principle, to characterize the replica 
symmetric region in terms of its own physical content, without resorting to 
features of the broken replica symmetry phenomenon. We will show that this 
is possible indeed. 

Our method is very simple. It can be considered a refinement of the 
interpolation method introduced in jB]. We will exploit also in a crucial 
way the established fact that overlap fiuctuations are Gaussian in the replica 
symmetric region, as shown for example in and |H]. 

The paper is organized as follows. In section 2 we recall the main fea- 
tures of our improved interpolation method. We consider an abstract general 
setting, and derive upper and lower bounds for the generalized free energy. 
These bounds are the essential ingredients to be exploited for the rigorous 
study of mean field spin glass models. In section 3 we introduce the mean 
field spin glass model in the replica symmetric region, and state the main re- 
sults. The quadratic replica coupling for two replicas is introduced in section 
4, and the main features of the corresponding replica symmetric region are 



investigated. Section 5 is devoted to the calculation of overlap fluctuations 
for systems containing two coupled replicas. These calculations follow the 
pattern outlined in jH] in the case of a single replica, and are a very simple 
extension of the treatment given in [Zj. The main results of this paper are 
proven in section 6, where it is shown in particular that the replica symmetric 
region extends up to the Almeida-Thouless line, as expected. We make also 
the comparison with the replica symmetric region given by Talagrand. In 
particular, we are able to isolate, in full generality, the parameter ruling the 
validity of a generic Ansatz for the free energy density. Section 7 is dedicated 
to some conclusion and outlook for future developments. 

2 Upper and lower bounds from interpola- 
tion methods 

For the sake of future reference, we give here a complete and self-contained 
treatment. Let us introduce the following general setting. Let t/, be a family 
of centered Gaussian random variables, i = 1, . . . ,K, with covariance matrix 
given by E(UiUj) = Sij. We treat the set of indices i now as configuration 
space for some statistical mechanics system, with partition function Z and 
quenched free energy given by 

Elog^Wiexp(Vtf/,) = E\ogZ, (1) 

i 

where Wi > are generic weights, and t is a parameter ruling the strength of 
the interaction. 

The following well known derivation formula plays a leading role in the 
interpolation method 

—E log '^ Wj expjVtUi) = 

i 

-E{Z'^ ^ Wi exp{ViUi)Sii 

i 

--E{Z-^ J2 Yl ^*^J' exp(v^t/*) exp{VtUj)Sij). (2) 

i 3 

The proof is straightforward. Firstly we perform directly the t derivative. 
Then, we notice that the random variables appear in expressions of the form 
E{UiF), were F are functions of the f/'s. These can be easily handled through 



the following integration by parts formula for generic Gaussian random vari- 
ables, analogous to the Wick theorem in quantum field theory, 

E{U,F)=Y,S.,E{^F). (3) 

Therefore, we see that always two derivatives are involved. The two terms 
in (J2I) come from the action of the Uj derivatives, the first acting on the 
Boltzmann factor, and giving rise to a Kronecker 5jj, the second acting on 
Z~^ , and giving rise to the minus sign and the duplication of variables. 

The derivation formula can be expressed in a more compact form by 
introducing replicas and suitable averages. In fact, let us introduce the state 
uo acting on functions F of i as follows 



uj{F{i)) = Z-^ywiexp{ViUi)F{i), (4) 



together with the associated product state fl acting on replicated configura- 
tion spaces ii, 12, ... , is- By performing also a global E average, finally we 
define the averages 

{F),^EniF), (5) 

where the subscript is introduced in order to recall the t-dependence of these 
averages. 

Then, the equation ^ can be written in a more compact form 

— E log"^ Wiexp{^/tUi) = -{Si,,,)t - -(^iiia)*- (6) 

i 

By following the same procedure, we can also easily establish the following 
general derivation formula for () averages 

"jT\-^(*l))t = '^{^hWiih ~ Si2i2 — 2Sij^i2 + '^^iiizllt- (7) 

Here, up to three replicas are involved. In fact, the first t derivative will let a 
second replica appear, while the subsequent d/dAj derivatives will introduce 
also a third replica. 

Let us now go to the interpolation and comparison arguments. Let Ui and 
Ui, for i = 1, . . . ,K, be independent families of centered Gaussian random 
variables, with covariance matrices given by E{UiUj) = Sij and E{UiUj) = 
Sij. For < t < 1, introduce the interpolating function 

a{t) = E\og^Wiexp{ViUi) exp{Vl - tUi), (8) 



so that we have the boundary terms 

a{0) = ElogJ2wiexp{U,) (9) 



a 



'l)=ElogJ2w^exp{U,). (10) 



By applying the t derivative as before, and taking into account that the 
1 — t term will bring an extra minus sign, we arrive to the basic interpolation 
formula 

— Q;(t) = -(oj^jj — Oi-^i^)t — -\bi^i2 — Oi-^i^Jt- (11) 

In the simplest case, the comparison argument runs as follows. Assume that 
the covariances satisfy the inequalities for generic configurations 

E{U,Uj) = S,, > E{U,U,) = S,j, (12) 

and the equalities along the diagonal 

E{U,U,) = Su = E{UA) = Su, (13) 

then we have 

|a(t) < 0, (14) 

with the obvious consequences 

a(l) < a(0), a{t) < a{0). (15) 

In the applications, a{l) is the quantity of interest, while a(0) is connected 
to some comparison system, and can be easily calculated. Therefore, we call 
p5|l an upper bound formula. 

Considerations of this kind are present in the mathematical literature of 
some years ago. Two typical references are Pj and ^^. 

If the stated conditions for the covariances ()12ll3p hold, then we also have 

|(a(0)-a(t)) = i(Ai2)t, (16) 

where 

Ai2 = S,,,, - 4,, > 0. (17) 

By integration, we have also 

1 r^ 

a(0) = a(l) + - / {Au)tdt. (18) 

^ Jo 
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We call "sum rule" an expression of this kind, by following the terminology 
introduced in jH]. 

The derivation formula and the comparison argument are not restricted 
to the Gaussian case. Generalizations in many directions are possible. For 
the diluted spin glass models and optimization problems we refer for example 
to jTT], and to [12j, and references quoted there. 

Let us show how to derive lower bounds, in the opposite direction. We 
follow the general strategy expounded in jHj, with some relevant improve- 
ment. 

Introduce a system made by two replicas with some appropriate coupling 
between them. Now the interpolating quenched free energy is defined as 

a(A, t) = 

-E log Y^ w,Wj exp{Vi{U, + Uj) + Vl - t{Ui + Uj)) exp{X{Sij - Sij)), 

(19) 

depending on two parameters, the interpolating parameter < t < 1, and 
the coupling parameter A > 0. We have also introduced the factor 1/2 in 
order to keep the same overall normalization. In fact, we have (5(0, t) = a{t), 
if A is put to zero. 

It is simple to calculate the t derivative by a direct application of the 
general formula (fTT|) . Of course, now we have to take into account that 
the original system is made of two (coupled) replicas, let us say (1,2). The 
replication process, therefore, will introduce a replicated system made by the 
replicas (3,4), coupled in the same way as (1,2). Now the diagonal term 
gives a nontrivial contribution, and we find 

(20) 





dta{X,t) = -{Ai2)x,t- 


- (Ai3)A,t, 


where as before 






Ai2 


= Si-^i^ — Si^i^ > 0, Ai3 


^ *-'nJ3 *-'«iJ3 — 0. 



(21) 

The terms A14, A23, A24 do not appear explicitly, because we have exploited 
the (1 — 2) and (3 — 4) symmetries. Of course, because of the A coupling, 
there is no symmetry of the type (1 — 3), and analogous. However, at A = 0, 
we have 

(Ai2)o,t = (Ai2)t = (Ai3)o,i. (22) 

On the other hand, the A derivative is given simply by 

9A«(A,t) = i(Ai2)A,t. (23) 



Introduce now a velocity field f (A, t), and trajectories t — > A(t), such that 
A(0) = Ao, and A(t) = v{X{t),t). Then we have for the total t derivative 

|a(A(t),t) = ^il + vim,t)){Au)x,t - (Ai3)a,. (24) 

In our previous paper P, we made the very simple choice v{X,t) = —1, 
so that X(t) = Xq — t. In this case 

^«(A(t),t) = -(Ai3)A,*<0. (25) 

at 

Therefore 

a(A(t),t) < a(Ao,0). (26) 

A more efficient way is to go along the level lines of a, with the natural 
choice 

viX,t) = -l + 2^^^, (27) 

SO that 

^«(A(t), t) = 0, «(A(t), t) = a(Ao, 0). (28) 

This choice will play a central role in our treatment. Notice that v{0, t) = 
1. Let us now derive the lower bound. From the definition ()19|) of a, by 
applying Jensen inequality, we easily get 

a{X{t),t) = a(t) + -Elogn(exp(A(t)Ai2)) 

> a{t) + ^{Au)t. (29) 

By recalling (fT^ and (PUJ) or (PHj). we obtain 

A(t) — (a(0) - a{t)) < 5(Ao, 0) - a{t). (30) 

Provided the condition A(t') > holds for < t' < t, we can easily integrate 
(IHUj). and derive the basic lower bound for a(t) 

0<a(0)-a(t))<(a(Ao,0)-«(0))(exp(y^ J^) - '^)- (31) 

In the following we will make extensive use of this lower bound. 



3 The mean field spin glass model in the replica 
symmetric region 

The generic configuration of the mean field spin glass model is defined through 
Ising spin variables ai = ±1, attached to each site i = 1,2, . . . ,N. 

The Hamiltonian of the model has a random character, and can be defined 
as follows. First of all let us introduce a family of centered Gaussian random 
variables IC{a), indexed by the configurations a, and characterized by the 
covariances 

E{}C{a)K:{a')) = q\a,a'), (32) 

where q{a, a') are the overlaps between two generic configurations, defined 
by 

g(a,a') = iV-i^a,a:, (33) 

i 

with the obvious bounds —1 < q{a, a') < 1, and the normalization q{a, a) = 
1. 

At a generic inverse temperature /?, for an external magnetic field of 
strength h, the (random) partition function is defined in the form 

Z^{P,h;IC)= J2 exp(/3y^/C(a))exp(/3/i^a,), (34) 

(Tl...cr]v i 

which will be the starting point of our treatment. 

By now standard methods [13j assure the existence of the infinite volume 
limit in the form 

lim N-^ E log Z N {(3, h;IC) = sup N~^E log ZN{f3,h;IC)=a{(3,h). (35) 

Now we introduce the Sherrington-Kirkpatrick replica symmetric solu- 
tion, by following the methods explained in [Hj. 

For a generic trial parameter q, and an interpolating < t < 1, define 
now 

a^it) = N~'E\og Y, exp(M^^/C(a))exp(/?^(/i + v/r^v^J,)a,) 

(T1...(TJV i 

+ l(3'{q-mi't), (36) 

where we have introduced independent unit Gaussian fields Jj, acting as 
magnetic fields on each site. Now E denotes average over the /C's and the 

J^'S. 
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For t = 1 we have 

aN{l) = N-^E\ogZM{(3,h-K,), (37) 

while for t = 0, the Boltzmannfaktor factorizes, the sum over a's can be 
performed exphcitly, and we obtain the rephca symmetric trial function 

a^(0)=log2+ I log cosh (3 {h + ./^z)dn{z) + ^p\q -If. (38) 

Here dfi{z) is the unit Gaussian measure. Of course, aiy{0) does not depend 
on A^, and will be written as a{0) in the following. 

By the same procedure as in Section 2, we can easily calculate the t 
derivative in the form 

laN{t) = -^{{qu-qf)t<0. (39) 

Clearly we have the uniform bounds 

ajv(t) < a(0), a^(l) < a(0). (40) 

In view of these bounds, since g is a free parameter, it is convenient to 
take into (jHH|) the value that minimizes a(0). This value q{f3,h) is uniquely 
defined, as shown in S\, and satisfies the implicit equation 

q{P,h) = [tanh\p{h+^/q{p,h)z) d^i{z). (41) 



It turns out that g(/?, /i) > 0, unless /i = and /5 < 1 (the annealed region), 
where g(/3,0) = 0. From now on, we assume that this minimizing value 
g(/3, h) has been chosen in all our formulae for q. In this case a(0) in (jHHj) 
gives the replica symmetric Sherrington-Kirkpatrick value asK{.l3, h) = a(0), 
and we have 

|(«5i.(/3, h) - a^it)) = ^{iq^2 - qf)t > 0. (42) 

Now the bounds in ()4()j) are given in the optimal form 

aN(t) < asK{P,h), aAr(l) < asK{P,h). (43) 

Let us now introduce the parameter yg = VoiP, h) defined by 

{(3,h)= [ cosh-\/3{h + ^/q{/3, h)z))dfi{z). (44) 



2/0 ( 

On the (/?, h) plane, P >0, h > 0, define the SK region, the AT line and the P 
region as the regions where P^yo{P, h) < 1, P'^i/q^P, h) = 1 and P'^yo{P, h) > 1, 
respectively. 

The main result of this paper is contained in the following theorem. 
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Theorem 1. In the replica symmetric region SK, where l3'^yo{l3,h) < 1, 
the Sherrington- Kirkpatrick value gives the right free energy in the infinite 
volume limit, i.e. 

lim N-^E log Zjv(/3, h; /C) = a(/3, h) = asxiP, h) = 

N—^oo 

f(log2+ flogcosh{(3{h + ^z))dfi{z) + ^(3\q-iy). (45) 

Let us now consider the interpolation function a^lt) as defined in (jHBj) . 
The methods explained in show that there exists t(/3, h), < t(/5, h) < 1, 
such that for any < t < i{(3, h), we have in the limit 

lim aAr(t) = asK{P,h). (46) 

We assume that t(/?, h) is the largest value with this property, and call the 
interval < t < i{f3, h) the replica symmetric interval for t. Notice that this 
interval is always nontrivial, since i{f3, h) is always larger than zero. 
The following Theorem gives a full characterization of i{(3, h). 

Theorem 2. In the replica symmetric region SK, where j3'^yo{l3,h) < 1, 
i{P,h) takes its maximum value i{P,h) = 1. In the broken replica region, 
where P'^yo{P, h) > 1, the t(/?, h) takes the value 

mh)=(3-'yoiP,h)-\ (47) 

Clearly, Theorem ^ follows from the first part of Theorem El because we 
can take t = 1 in (jlH|l and recall (jSZj)- 

On the other hand, the first part of Theorem |21 follows from the sec- 
ond part, by taking the limit P'^yo —>■ 1. The argument by Toninelli [Ti] , 
based on the broken replica symmetry bounds of |^, shows that, in the bro- 
ken replica symmetry region, the parameter i{f3, h) can not be greater than 
(3~'^yo{f3, h)~^ . Therefore, in order to prove TheoremEl it is enough to show 
that i{(3, h) can not be less than (3~'^yQ{(3, h)~^. This will be our objective in 
the following. 

4 Quadratic coupling of replicas 

Now we couple two replicated systems, with spin variables a^^\ a^'^\ in the 
same spirit as in jH] and in Section 2. In analogy with ()19|) the auxiliary 
function a^i^, t) is now defined, for A > and < t < 1, as 

«iv(A,t) = ^Elog Y, exp(fii)exp(52)exp(^A/?2iV(gi2-g(/?,/i))') 

+l(3\q-ini-t), (48) 
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where Bi and B2 are the two rephcas of the quantities appearing in the 
Boltzmannfaktor of aj\[{t) in (jHBj) . 

The function aiv(A,t) is convex increasing in A for each value of t, and 
we have 

ajv(0, t) = aAr(t), Q;Ar(A, t) > aAr(t). (49) 

The t and A partial derivatives are easily calculated in analogy with (|^n|l 
and ^ 

d,a^{\, t) = ^((gi2 - qr)x,t - y ((gi3 - qf)x,t, (50) 

d,a^{X,t) = f^{{qu-q)')x,t. (51) 

Let us go to the infinite volume limit A^ -^ 00, by exploiting for example 
the methods outlined in p[S], and call a(A,t) the resulting limit of aAr(A,t). 
For any < t < t(/3, h), let A(t) be the highest value of A for which a{X, t) = 
C(SK{(3,h). We call the region < t < i{(3,h), < A < X(t) the rephca 
symmetric SK region in the (A, t) plane. 

In 10], the value of A(0) was calculated explicitly through a variational 
principle. In particular, in the case h = 0, q = 0,it turns out that A(0) = /9~^. 
Moreover, it was shown that the triangular, or trapezoidal, region < t < 
min(l,A(0)), < A < A(0) — t, is contained into the replica symmetric 
region. Another interesting property of the replica symmetric region is the 
following. Assume that (Ai,ti) is in the SK region, then all points ti < t < 
min(l, ti + Ai), < A < Ai — t + ti, are in the SK region. All these properties 
follow, according to [H], from the bound (jSH), adapted to the aN{t) defined 
in ()36|1 . with the simple choice v{X,t) = —1, by taking the infinite volume 
limit. 

A very important consequence is the following 

Theorem 3. In the case where t{(3,h) < 1, the following hound necessarily 
holds for X{t) 

X{t) < t{p, h) - t, (52) 

for anyO<t< i{P,h). 

Therefore, if t{f3, h) < 1, the replica symmetric region must shrink to the 
single point A(t(/3, h)) = at the end point t(/5, h). The proof is straightfor- 
ward. Should (jn21) fail, there would be a point ti where A(ti) > t(/?, h) — ti. 
But in this case, by the general properties pointed before, it would be pos- 
sible to enlarge the replica symmetric region beyond t(/?, /i), which is clearly 
impossible due to the very definition of t(/?, h). 

There are two important properties of the replica symmetric region that 
we immediately point out. 
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Introduce, for generic replicas s, s', the rescaled overlaps defined as the 
random variables 

U' = ^{qss'-q{(3,h)), (53) 

distributed according to the (. . .)A,t state. Then we have 

Theorem 4. For any (A,t) in the replica symmetric region, in the infinite 
volume limit, the ^ become centered Gaussian random variables, whose covari- 
ances can be explicitely calculated. These covariances may become singular 
only along the boundary line X, or at the boundary point i{f3, h), according to 
the various situations analyzed in the following. 

The proof is very simple. In fact, it is a simple adaptation of the same 
procedure exploited in |7j, and 0, for the case where A = 0. The effective 
calculation of the covariances exploits methods introduced in |H], [7j, and 0. 
We sketch the extension to the case where A > in the next section 5. 

In the replica symmetric region, we have not only the convergence of the 
a;Ar(A,t) auxihary function to asxiP, h), but also a stronger statement. 

Theorem 5. For any (A,t) in the replica symmetric region, with the possible 
exclusion of the end point i{f3, h), the infinite volume limit 

lim N{aNiX, t) - asKi/3, h)) = 0(A, t) (54) 

is well defined and can be explicitly given in terms of integrals of fluctuations 
as follows 

p' r ...... , p' '" 



0(A,t) = -^/ A{t')dt' + '-l A{X',t)dX', (55) 

where 



4 ./o ^ ' ^ JO 



A{t') ^ lim {^uh, A{X', t) ^ lim {e^2)x',t. (56) 

The function (j) could take the value oo at the boundary line X, according to 
the various situations analyzed in the following. 

The proof is very simple. In fact, we recall that ajv(0,t) = aN{t), and 
write 

N{aN{X, t) - asKiP, h)) = -N{asKiP, h) - aN{t)) + A^(«7v(A, t) - «;v(0, t)). 

(57) 
Then we exploit (jl^ and (f^ . By integration we get 

N{asK{P.h)-aM{t)) = L j (^^f^)^,dt', (58) 

iV(a^(A,t)-a,v(0,t)) = y /" (^12)^,*^^. (59) 

12 



By taking the infinite volume limit we get ()55|) . and Theorem |31 is proven. 

It is also convenient to introduce the level lines of the function 0. In fact, 
we follow the procedure exploited in section 2, by starting with the level lines 
for the the difference aAr(A,t) — asK{[3, h)). It is immediately recognized, as 
in section 2, that, for finite A^, the corresponding velocity field is given by 

v{K t) = -1 + 2i^^l^^i£^ ^ -1 + 2|M^. (60) 

((gi2 - qr)x,t (Cl2)A,t 

In the replica symmetric region, fluctuations can be completely controlled. 
Therefore, in the infinite volume limit, we have that the velocity field of the 
level lines for the limiting function are given explicitly by 

where A{\,t) is defined in fl^ . and 

DiX,t)= ]\m{e^s)x,t■ (62) 

Notice that A{X, t) and D{X, t) are different in general, because the A coupling 
destroys the (2 — 3) symmetry. They are explicitly calculated by following 
the method outlined in section 5. Of course, at A = 0, we have A(0,t) = 
Z}(0,t) = A{t). In particular f(0,t) = 1, for any t < i{P,h), where A{t) is 
surely not singular. 

Let us give a simple example. We consider the annealed case, where 
h = 0, q = 0. Here everything can be explicitly calculated. In particular, the 
methods of ^ give i{p) = 1 if /3 < 1, and i{p) = /^-^ if /3 > 1. The annealed 
region is given by < t < t(/5), < A < A(t), where \(t) = /?"^ - 1. In the 
annealed region, by following the methods outlined in section 5, we find that 
the variances are given by 

A{t) = l/{l-tf3''), D{X,t)=A{t), A(A,t) = l/(l-(A + t)/52). (63) 

Notice that A{X,t) becomes infinite on the A line, where v{X{t),t) = —1. In 
general we have for the velocity field, from (J6T| and (J63|l . 

The level lines are immediately found by integration of the streaming equa- 
tion A(t) = v{X{t),t). Let Ao be a parameter with the condition Aq < Aq = 
/3~^. Then the family of level lines is given by 

A(t) = ir' -t)- ifr' - Ao)(l - Phf, (65) 
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where in any case it must be < t < 1 and A(t) > 0. If Xq > 0, then the 
trajectories start from Aq at t = 0. If Aq < then the trajectories start from 
A = at t = Ao/(l — /S^Ao). The value of the function 0, as defined in (jUj), 
can be easily found from (jHH|l in the form 

where Aq is the parameter characterizing the trajectory passing through (A, t). 

In this simple completely solvable annealed case, we can point out the key 
phenomenon appearing on the (A,t) plane. Assume (3 > 1. Then we know 
i{l3) = l3~'^ < 1. The annealed region is the triangle < t < t(/3), < A < 
/?~^ — t. At the end point t = t(/5), A = 0, the variance A(t) becomes infinite, 
the function 0(O,t(/9)) is not well defined, and moreover all trajectories have 
there an osculating focus with common slope f (0,t(/3)) = —1. 

We will show in the following that this kind of picture holds qualitatively 
also in the general replica symmetric case. 

5 Calculation of overlap fluctuations 

Overlap fluctuations have been calculated in the infinite volume limit in [S] 
and jZj, see also jH]. 

Firstly, we consider the case A = 0. The initial conditions at t = follow 
easily from the central limit theorem. In fact we have 

Bo = lim (6263)0 = q-f, (67) 





Ao -= 


E lim(e?2)o = l- 

N—*oo 


-f, 




Co -= 


= lim (6264)0 = 


/ ta 


is 


also convenient to define 








yo = Ao- 2Bo 


+ Co 



tanh\p{h + ^/q{p,h)z))dfx{z) - f. {6t 



Zo = Ao- ABo + 3Co. (69) 



It is immediately recognized that yo has the value already given in 
Moreover, the bound Zq < j/o holds. 

In the replica symmetric interval, in the infinite volume limit, the ^ be- 
come centered Gaussian random variables, with explicitly given covariances 
(see jH], [Z!, and [Sj) according to 

A{t)= limjv-.oo(6'2)t =?>y{t)~2z{t)+w{t), (70) 

B{t)= lim^_^oo(6263)t =^y{t)-^z{t)+w{t), (71) 

C{t)= lim,v-.oo(6264)t =y{t)-z{t)+w{t), (72) 
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where 

In order to calculate the covariances for A > we follow a very simple 
strategy. For any t in the replica symmetric interval, we take A{t), B{t), C{t) 
as initial conditions, and consider the streaming equations in A. These are 
easily solved in the infinite volume limit. There is a mild complication. Since 
the symmetry between replicas is explicitly broken by the A couplings, there 
are more basic covariances to be taken into account. Let us consider some 
examples. 

Define, in agreement with (|HF)|l . 

A(A, t) ^ hm {e,2)x,t, C(A, t) ^ lim {^,2U)x,t, (74) 

N^oo Af— >oo 

with the initial conditions A{0,t) = A{t), C{0,t) = C{t). At fixed t, the 
streaming equations are 

J^(e?2)A,* = Y^ei2ieu-eu))Kt, (75) 

d l3^ 

— (62^34) A,t = -^(6264(^12 + ^34 ~ ^56 ~ ^78))A,i- (76) 

Now we go to the infinite volume limit, take into account the Gaussian na- 
ture of the resulting fluctuations, and exploit the partial symmetry among 
replicas. We find the equations 

A'{X, t) = A{X, tf - C(A, tf, C\\, t) = 2A{X, t)C{X, t) - 2C(A, t^, (77) 

where the prime denotes the partial derivative with respect to the rescaled 
parameter A/3^. These equations, with the given initial conditions, are easily 
integrated in the form 

^(^•*) = 11 - xp^A{t) - C{tW '"> 

Analogous procedure can be followed for the other covariances. For example, 
for 

B{X,t) = lim {^u^i3)x,t, D{X,t) = lim (^2^)^,*, ^(A,t) = lim {^i3^,e)x,t, 

N^oo N—*QO N^oo 

(80) 
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we find the equations 

B' = AB + BC-2CE, D' = 2B^-2E'^, E' = 2BC + AE ~3CE, (81) 

which are easily solved by taking into account the initial conditions -6(0, t) = 
Bit), D{0,t) = A{t), E{0,t) = C{t). 

We end this section with the following qualitative estimates, which are 
useful for the investigation of the properties of the level lines introduced in 
the previous section. 

Theorem 6. The following estimates hold 

^(A,t)<^-^^^, D{X,t)>A{t), v{X,t) >1-2X(3' A{t). (82) 

The proof is simple. By neglecting the C^ term in the first equation 
in (fTTjl . we obtain A' < A"^, which gives the first inequality in (jH^ upon 
integration. On the other hand, by exploiting ()81|1 . we find 

(B-E)' = {B-E){A-C). (83) 

In general A > C hj Schwarz inequality, and B > E a.t X = 0, because here 
it reduces to the obvious B{t) > C{t). Therefore B > E for any A. But 
then, from ()81|) . we derive D' > 0, and establish D{X,t) > A{t). Finally, the 
lower bound for v comes from the previous bounds, and the definition ()6ip . 

6 The replica symmetric region extends up 
to the Almeida-Thouless hne 

Now we are ready to prove Theorem El Assume that P'^yo{P,h) > 1. We 
show that i{P,h) must have exactly the value i{P,h) = P~^yQ{(3,h)~^. By 
Toninelli argument in [T3], we already know i{/3, h) < j3~'^yQ{j3, h)~^. Suppose 
i{(3, h) < /?~^|/o(/?, h)~^. Then the variance A{t) would stay finite up to the 
point t = i{P, h) included, as the explicit solution in (fTOI) shows (recall that 
Zf) < yo)- As a consequence, the function 0, defined in fj5l|l . would be well 
defined at the point (A = 0,t = i{P,h)), due to the integral representation 
()55|) taken at A = 0. Here we find an immediate contradiction. In fact, 
due to the condition (fS^ . in the point (A = 0,t = i{P,h)) all level lines 
should converge. But each of them carries its own value for the function 
0, depending on the initial condition Aq. For example, along the trajectory 
starting from A = at t = 0, we must have = 0. As a general fact, the 
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level lines can not focus on a point where is well defined. We conclude that 
i{P, h) < l3~'^yo{l3, h)~^ is impossible. We must have t(/3, h) = (3~'^yQ{l3, h)~^. 
Therefore, Theorem |2l is proven, together will all results of this paper. 

We end this section with an argument showing the unwillingness of the 
level lines to focus in a point where A{t) is not singular. Let us recall that 
we have proven the bound f (A,t) > 1 — 2X13"^ A{t) (see (jH^ ). for the velocity 
field of the level lines. Now assume that A{t) is finite along the whole replica 
symmetric interval, so that A{t) < a, for some positive constant a. Then we 
would have v(A, t) > 1 — 2A/5^a. This can be easily integrated. For example, 
we would find that the level line starting from A = at t = 0, should stay 
above the line 

A(t) = (1 - exp(-2t/52a))/2/32a, (84) 

which is clearly incompatible with ()52|) . 

In conclusion, our general argument shows that, on the (/3, h) plane, the 
boundary between the replica symmetric region and the broken replica sym- 
metry region is characterized by a singularity in the rescaled overlap fluctua- 
tions, exactly as happens in the well known case of vanishing magnetic field. 
Therefore the transition is exactly at the Almeida-Thouless line. 

It is interesting to compare our results with the characterization of the 
transition line given by the very clever method of Talagrand |3j. Talagrand 
exploits broken replica symmetry bounds for a system made by two coupled 
replicas. As a result, the transition line is given in terms of the onset of 
the instability of the replica symmetric solution in the frame of a one step 
breaking of the replica symmetry in the Parisi representation. Therefore, 
the problem of comparing the Talagrand transition line with the Almeida- 
Thouless line is open. In our method, we do not employ any information 
coming from the breaking of the replica symmetry, but exploit only the pe- 
culiar properties of the quadratic replica coupling. It is also important to 
remark that Talagrand method has been extended to give the proof of the 
validity of the Parisi representation for the free energy in the infinite volume 
limit 0. In next section we deal with the problem of extending our method 
to more general situations. 

7 Conclusion and outlook for future develop- 
ments 

We have seen that the control of the rescaled overlap fluctuations and the 
peculiar properties of the quadratic replica coupling allow to show that replica 
symmetry holds up to the Almeida-Thouless line. No information coming 
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from the pattern of the spontaneous rephca symmetry breaking phenomenon 
is necessary. 

In its very essence, this method allows to isolate an important parameter 
ruling the phenomenon. The replica symmetric representation holds, pro- 
vided the parameter P'^{Ao — 2Bq + Co) is less or equal to 1. Therefore, the 
ruling parameter is given in terms of fluctuations at zero coupling! 

It is important to remark that the flow equations for fluctuations do not 
depend on the particularly chosen representation. The solutions do depend 
on the representation, only through the initial conditions {Aq,Bo,Co). 

Therefore, it is possible to extend our method also to other represen- 
tations. In particular, in the frame of the Aizenman-Sims-Starr extended 
variational principle P£j, we can take the representation based on Derrida- 
Ruelle-Parisi probability cascades (see for example [H]). Here sum rules also 
hold, and we can consider a generalized quadratic replica coupling between 
replicas. We plan to report on these problems in a future publication. 
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